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Abstract 

In this paper, starting from the microscopic dynamics of isolated dislocations, we explain how 
to derive formally mean field models for the dynamics of dislocation densities. Essentially 
these models are tranport equations, coupled with the equations of elasticity. Rigorous results 
of existence of solutions are presented for some of these models and the main ideas of the proofs 
are given. 
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Introduction 

In this paper we are interested in several mesoscopic models involving the dynamics 
of dislocations. These models are important for the understanding of the elasto-visco- 
plasticity behaviour of materials. For crystals, at the microscopic level, the origin of 
plasticity is mainly due to the existence of defects called dislocations. Dislocations are 
line defects that can move in the crystal when a shear stress is applied. On the other 
hand each dislocation also creates a stress field. This leads to a complicated dynamics 
of these defects that we consider below. See Figure Q] for an example of complicated 
pattern of dislocations in real crystals. 

Organization of the paper. 

Section [Q presents the microscopic modelling for the dynamics of dislocation straight 
lines. In Section [2j we explain how to derive formally a two-dimensional mean field 
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Figure 1: Dislocation lines in Nickel Aluminium alloy. 



model, called the Groma, Balogh model. A rigorous result of existence of solutions is 
given and other results are also given for a one-dimenional submodel. In Section [3l the 
Groma, Czikor, Zaiser model is studied. In particular is presented a simulation of the 
deformation of a slab under an external shear stress. In Section the formal derivation 
of a mean field model for densities of dislocation curves is presented. Concluding remarks 
are given in Section [5l 

1 The microscopic modelling 

1.1 Preliminary: the stress field created by an edge dislocation 

In the space coordinates x = (xi,x 2 ,x 3 ) with basis (ei,e 2 ,e 3 ), we consider the case a 
material filling the whole space and with a single dislocation line which is the axis x 3 . 
To this line is associated an invariant, called the Burgers vector b. In the case of linear 
isotropic elasticity with Lame coefficients A, /i (satisfying /i > 0, 3A + 2/i > 0), the stress 
created by the dislocation is given by 



where H is the Heaviside function, 5 is the Dirac mass, n = e 2 is the normal to the 
slip plane and (mi,m 2 ,m 3 ) is the three-dimensional displacement field. In the special 
case where b = e\, the dislocation is called an edge dislocation (the Burgers vector is 
perpendicular to the dislocation line). In this case we have 




(1.1) 



where 





(1.2) 
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The stress is assumed to satisfy the elasticity equation of equilibrium 

0, j = 1,2,3. 



i=l,2,3 



dxi 



(1.3) 



The solution to this equation is known since the work of Volterra [16j. We have (see 
Hirth, Lothe |7|): 

/ x 2 (3x 2 + x 2 ) Xi(x 2 — x 2 ) \ 



(Tlx) 



[x\ + x|) 2 



(xf + x 2 ,) 2 







1 1 *^ 2 ) ^ 2 1 ^ 2 ) 



(xf + x 2 ) 2 



(xf + x 2 ,) 2 



with 



2tt(1 - u) 



> where the Poisson ratio is v 



0/ 



A 



2(A + /i) 



1 

4 '2 



:i-4) 



1.2 A two-dimensional microscopic model 

In this section we describe very formally a microscopic model describing the dynamics 
of dislocation lines in the particular geometry where all lines are parallel to the axis X3. 

Considering the cross section of these lines, we can reduce the problem to a two- 
dimensional problem where each dislocation line can be identified to its position (xi, x 2 ). 
For % — 1, N, let us call X % the positions in M 2 of these dislocations (see Figure [2]). 
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Figure 2: The cross section of the dislocation lines. 

We consider moreover the particular case where each dislocation can only move along 
the direction e%. Those dislocations X 1 are known to be characterized by an invariant 
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called the Burgers vector b{ = Ei&\ where 

Si = ±1. 

Finally the dynamics of these dislocations is 

dX\t) 



dt 



h a 12 (X l (t)) } for i=l,...,JV. (1.5) 



Here a 12 is one component of the stress. In the case of linear isotropic elasticity, it is 
known that 



crii{x)= ^ £ i a o(x-X l ), where a (x) = a ( :2y2 . 



i=l,...,JV 



where <To is the (1 — 2)-component of the stress given in (11.41) . Here ao(x) appears to be 
the shear stress created at the point x by a single dislocation positionned at the origin, 
with Burgers vector t\. With the convention that the self-stress created by a dislocation 
on itself is zero, i.e. formally cr (0) = 0, we see that we can rewrite the full dynamics of 
particles satisfying X % 7^ X 7 for % ^ j, as follows 



dX l 

-dT = ei 



e t eM xt - xj ) ) : fOT i = l,...,N. 

y je{i,...,N}\{i} 



Let us now introduce the "density" of dislocations 

d ± (x,t)= Yl So(x-X\t)) 

ie{l,...,N} : Si=±l 

where 9 + and 6~ correspond respectively to the dislocations with positive and negative 
Burgers vector. Again with the convention cr (0) = 0, we can rewrite the dynamics as 

df + div (± J 0±) = with J = ei a *(6 + -6-). (1.6) 



2 Mesoscopic models 

2.1 Formal derivation of the two-dimensional mesoscopic model 

A natural question is what is the effective model at the mesocsopic scale corresponding 
to the microscopic model (11.61) ? It seems to be an open and difficult question. As a naive 
answer, we can consider as a first candidate to the mesoscopic model, the mean field 
model corresponding to the model (11.61) . where now each density 6 + and 9~ can be seen 
as a "continuous" quantity. In this mean field model, we have in particular neglected 
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the short range dynamics. One convenient way to rewrite the model is to introduce one 
primitive p of the densities 9 ± such that 

p±(x,t) = 6 ± (x,t)>0. 

Then we can rewrite (11.61) as (with a zero constant of integration): 



pf = = F cr 12 P% 



(2.7) 



with 



ffi2 = (go)*i * (P + - P )• (2-8) 

System fj2.7D - (j2.8p is an integrated form of the Groma, Balogh model (11.6P (see [5]). It 
turns out that by construction a solves: 



i=l,2,3 

(Tij = 2peij + A I ^2 e kk\ 

\fc=l,2,3 / 



'i.l- 



(2.9) 



1 ( dUj dUi i , n 



2 V dxi dxj 



where e° is given in (11.21) . Here (12.81) appears to be a representation formula for the 
solution of the equations of elasticity (12.91) . 



2.2 Existence of a solution in the two-dimensional periodic case 

In the case where the dislocations densities 6 ± {x, t) are Z 2 -periodic in x G M 2 such that 

dx\ 9 (xi, X2,t) = L > 0, 



it turns out that the representation formula (12.81) can be rewritten 

cr 12 = a RlRl(p + — p~) with a = 4/i(A + p)/(\ + 2p), 



(2.10) 



where Ri for 2 = 1,2 are the Riesz transform given in Fourier series coefficients on 

T 2 = M 2 /Z 2 by 



Ck(Riif)) 



^-c k (f) for all A;=(A; 1 ,A; 2 )GZ 2 \{(0,0)} 







where 



Ck(f) 



if k = (0, 0) 

e~ 2ink - x f(x) d 2 x. 



T 2 



We will also make use of the norm on the following Zygmund space 
||/||zW)= inf| 7 >0, Jj£]n(e + ^\<1 

We assume that 

P ± (',0) = p ± 

and we make the following assumptions on the initial data 

Pq {x\ + m, %2 + I) = mL + p Q (xi, x 2 ) for any m,l 

< (Pok>0 

k ||(Po)^'ilUlnL(T2) < C 

where C > is any fixed constant. 
Then we have: 

Theorem 2.1 (Global existence of a solution in the periodic case, [2]) 

Under assumption $2.12\) there exists a function (p + ,p~) which is a solution of \2. 7\) . 
$2.101) in the sense of distributions and with initial condition $2.12\) . such that p ± (-,t) 
satisfies $2.12\) for all time. Moreover p ± G C([0, +oo); L^IR 2 )) and 

a 12 eLl c ([0,+oo);Hl c (R 2 ))- (2-13) 

Remark 2.2 Recall that there exists a constant C > such that the general inequality 
holds for functions f G if^T 2 ) and g G LlnL(T 2 ) 

II/^HlHT 2 ) < C||/lk 1 (T 2 )llfi , IUlnL(T2)- 

Then the product o~\ 2 ■ {p ± ) Xl in $2. 7\) is defined because of both estimate $2.13\) on cr 12 
and estimate L°°([0, +oo); L lnL(T 2 )) on (/7 ± ) X1 - 

Indeed, in the proof of Theorem 12.11 the main tool is to consider the entropy 
5 W = E/ ± (;t)\n6 ± (-,t) with 9 ± = (p ± ) Xl 

which satisfies formally the following a priori estimate 

S(t) + a [ [ (R 1 R 2 (e + - e~)) 2 < S(0). 

Jo Jt 2 

The uniqueness of the solution remains an open problem. 
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) 

(2.11) 

G Z 

(2.12) 




Figure 3: A 1-D sub-model invariant by translation in the (—1, 1) direction. 



2.3 A one-dimensional mesoscopic submodel 

In this section, we present a submodel where the uniqueness of the solution is known. 



Let us now consider a solution (p + , p~) of (12.71) . (12.101 ) only depending on the variable 
y = X\ + x 2 and on the time t (see Figure [3]). 

In that case, we can rewrite the system of equations for (p + (y,t), p~(y,t)) as follows: 



Pi 



Pt 



-Ci|(p H 

cA(p + - 



p ) + c 2 / dz (p + (z, t) - p (z, t)) } p 



p ) + c 2 J dz (p + (z,t) - p (z,t))\ p y 



(2.14) 



with the constants 



Cl 



p(\ + fi) 



> 0, c 2 



> 0. 



X + 2p X + fi 

For p^ > 0, this system looks very much like the Burgers equation in the rarefaction 
case (where no shocks are created). For this system, the following uniqueness result is 
available: 

Theorem 2.3 (Uniqueness for the one-dimensional submodel, [4]) 

Assume that the initial data p$ is Lipschitz non- decreasing and that for some L > 0, the 
functions y i— > p^{y) — Ly are 1 — periodic. Then there exists a unique viscosity solution 
(p + ,p~) to the system (i2.14\ ), f2.11\) . Moreover this solution is globally Lipschitz in space 
and time. 

In this theorem the notion of viscosity solution is the one introduced by Ishii, Koike 
[T2] for systems which are quasi-monotone. Roughly speaking, this is related to the fact 
that for c 2 = 0, this system has a comparison principle. More precisely, if p^ 1 and p ± ' 2 
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are two solutions of system (12 . 141) with C2 = such that 

P + 'W) < P + '\y,t) and p-'\y,t) < p-' 2 {y,t) for all y eR 

holds at time t — 0, then this is true for all time t > 0. 

Moreover for system (12 . 14H . it is possible to write an upwind scheme and to prove a 
Crandall-Lions type discrete-continuous error estimate for this scheme, as it is done in 
El Hajj, Forcadel [4]. 

Let us mention that an existence and uniqueness result for system ( 12.141) has also 
been obtained by El Hajj [3] in the framework of Hl oc (M) initial data with solutions in 

fl 1 to (Rx[0,+oo)). 

This system has also been studied in the case of periodic external applied stress. In 
the system (12.141) . this corresponds to add a time-periodic term to the quantity p~). 
Then for this non-local system, it is shown formally in Briani, Cardaliaguet, Monneau [I] 
(see also Souganidis, Monneau [15] for a local system in the stationary ergodic setting) 
that the long time behaviour of the system is an equivalent quasilinear diffusion equation. 



3 The GCZ one-dimensional mesoscopic model 



3.1 Presentation of the model 

In this section we consider the case where the three-dimensional material is a slab fl 
(-1,1) xl 2 with boundary conditions 



a ■ n 



±re2 for 

where r 6 1 is a fixed constant (see Figure [3J . 

e 2 



Xi = ±1 



(3.15) 
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Figure 4: Geometry of the slab. 

In that case, we can check that the solution to the equation (12.91) on Q supplemented 
with boundary conditions (13.151) . is 

ai2 = r and cr^- = if {i, j} ^ {1, 2} . 
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And then from the evolution equation (12.7ft . we see that in the case r > 0, the positive 
dislocations move to the right and the negative dislocations move to the left. In order 
to take into account the short range dynamics with accumulations of dislocations on 
the boundary of the material, Groma, Czikor and Zaiser [6] have proposed to modify 
equation (|2.7I) into the following equation (GCZ model) 

pf = =F5-i2Px! with (7i2 = °~i2 + n (3.16) 

where r& is the back stress created by the concentration of dislocations 

(6» + - 6~) X1 ± ± 

n = -Do e+ + e - Wlth 9 = P X1 

where D is a diffusion coefficient that we take equal to 1 to simplify the presentation. 
Introducing the quantities 

p = p + — p~ and re = p + + p~, 
the full GCZ system can be rewritten with y = X\ E I = (—1, 1) as 



Pt = Pyy ~ TK y 
K t K y = PtPy 



on Jx(0,+oo) (3.17) 



with boundary conditions on dl: 

A nfl +\ — n 

for all t > 0, (3.18) 



p(-l,t) = Q and p(l,t) = 
K,(—l,t) = —c and n(l,t) = c > 

for some constant Co and the initial conditions 



p(;0)=p° 

k{-,0) = re . 



(3.19) 



The non-negativity of the densities 6 ± at the inital time is equivalent to the following 
condition 

n° y > \p° y \ on /. (3.20) 
Moreover, we will assume the following condition 

(Py, Ky) is in C°°(I) and with compact support in /. (3-21) 



3.2 Main results 

For this system, we have the following result 

Theorem 3.1 (Global existence of a solution; jllj ) 

Assume that the initial data (p°,K°) satisfies $3.20\ )- [3.21\) . Then there exists two func- 
tions p E C 1 (I x [0, +oo)) and re E C° (7 x [0, +oo)) solution of ^l^ - ^M - ^TM) . 
Moreover for all time t>0, {13.20$ is satisfied with (p°, re ) replaced by (/?(•, t), re(-, t)). 
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Remark that in this theorem the notion of solution to (13.171) is the following. The 
first equation of ( 13.171) is satisfied in the sense of distributions, while the second equation 
of ( I3.17p satisfied by k is satisfied in the viscosity sense. This makes sense because the 
right hand side p t p y of the equation is continuous, and because the time derivation n t is 
multiplied by the non-negative quantity n y . 

Here condition (13.211) appears to be a technical condition to get the result. Althought 
this one-dimensional system ( 13.171) seems very simple, its mathematical study and the 
proof of Theorem 13.11 is particularly difficult. 

In the case r = 0, it is possible to get a uniqueness result. 

Theorem 3.2 (Uniqueness of the solution in the case r = 0; |9j) 

In the case r = 0, let us assume that the initial data (p°, k°) is Lipschitz and satisfies 
for some 5 > 0: 

«J > ^ 2 + (P° y ) 2 on I. 
Then there exists a solution (p,K) to $3. 11\) - [3.18\ )- [3~1~&) satisfying 



Ky(-,t)> ^S 2 + p 2 y (.,t) for all t>0. (3.22) 

Moreover this solution is unique among the solutions satisfying ( fg. 22\) . 
Main idea for the proof of Theorem 13.21 

The proof of Theorem 13.21 is based on the fact that ^J5 2 + p 2 is an entropy subsolution 
of the conservation law satisfied by v = K y , namely 

f n(y,t) \ 
v t = with n = ptPy 



Then the comparison principle between entropy solutions and entropy subsolutions shows 
that *),«(•,*)) satisfies (13T22D for all time t > 0. 

Main ideas for the proof of Theorem 13.11 

In particular, a first try shows that 

M= Ky~ \p y \ 

satisfies formally 

M t = a x M y + a M 

with 

a x = rsgn(p y ) - and a - ^ ^ Sgn ^^ 



(k y ) 2 ( K y) 2 K y 

By a maximum principle argument, we see that in order to guarantee that M > is 
true for every time, we have somehow to control the L°°-norm of p yyy /K y . This seems 
hopeless, because we would need to control moreover K y > from below. The idea is 
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then to replace system (13.171) by a suitable regularized system, which is the following for 

e > 0: 

Pt = (l+e)Pyy - TKy 

* = SKyy + P -^-r Py (3-23) 

Ky 

and the result is then obtained, passing to the limit e — > 0. Even the regularized system 
(13.231) is still difficult to study, because of the division by K y (see [TO]). But for the 
regularized system (13.231 ). it is possible to replace M by 



M1 = K y (;t)-y/jZ(t) + (p y (;t)y 

and to prove that M 7 > while the non-increasing function ^(t) satisfies 

^<-C E {l + \\ Pyyy (;t)\\L^I)) (3.24) 

where the constant C e blows-up as e — > 0. The striking remark is that to show in the 
case e > that M 7 is non-negative, we only need a control on the L°°-norm of p yyy , 
while in the case e = 0, in order to show that M is non-negative it was necessary to 
control the L°°-norm of p yyy /^ y which was worse (and not sufficient to conclude). 

On the other hand, proving a parabolic version of the well-known Kozono-Taniuchi 
inequality (see [13]) and using heavily the regularity theory for parabolic equations, we 
can prove for the regularized system (j3.23K that 

\\Pyyy(;t)\\L^l) < Ce Ct (l + | ln 7 (t) |) (3.25) 

for some constant C > also depending on e, among other things. Putting estimates 
(13.241) and ( 13.250 together, we see that we get the three-exponential estimate 



K y (-,t) > 7(t) > e 



-e 



for some constant c > depending in particular on e. Then the global existence of a 
solution to the regularized system (13.231) follows. Finally we recover a solution to the 
original system (13.171) taking the limit e — > 0. 

3.3 Numerical simulations 

Comming back to the system of elasticity, it is possible to compute the displacement. 
We find Ui = u 3 = and 

t r y 

U2{y,t) = -y+ dz p(z,t) 
P Jo 

In Figure [5], we show successively the initial state of the crystal at time t — without any 
applied stress, then the instantaneous (elastic) deformation of the crystal when we apply 
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the shear stress r > at time t = + . The deformation of the crystal evolves in time and 
finally converges numerically to some particular deformation which is shown on the last 
picture after a very long time. This kind of behaviour is called elasto-visco-plasticity 
in mechanics. Moreover, on the last picture, we observe the presence of boundary layer 
deformations. This effect is directly related to the introduction of the back stress r& in 
the model. 




a)t = 0~ b) t = 0+ c) t = +00 

Figure 5: Deformation of a slab for model ( 13.231) 



4 General dynamics of curved dislocations 

4.1 Preliminary on the dynamics of dislocation curves 

At sufficiently low temperature, dislocation curves are contained in the crystallographic 
planes of the three-dimensional crystal and can only move in those planes. Let us 
consider a closed and smooth curve T t moving in the plane (2/1,2/2)- The evolution of 
this dislocation can be modelled by a dynamics with normal velocity c(y, t) (see Figure 

ED. 




Figure 6: Schematic evolution of a dislocation line Tt by normal velocity c between the 
times t and t + At. 
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4.2 Transport formulation of the dynamics 

If we want to describe the dynamics of a high number of dislocations curves, it is inter- 
esting to describe this dynamics with a single quantity and a single equation. This is 
the goal of this subsection that we describe heuristically. 

Let us consider a smooth closed curve r in the plane (yi, 1/2) with basis (ei, 62). The 
curve To can be parametrized by its curvilinear abscissa s, such that 

To = \J{y(s)}. 

s 

Let us define for 6 G M/(2ttZ) 

n{6) = (cos#)ei + (sm9)e 2 and t{6) = (sin6 l )ei — (cos6 l )e2. 
We also define the angle 9(s) of the tangent to T at the point y(s) by 

|(.)=r W .)) 
and the curvature K^ (y) of r at the point y by 

Kr (y(s)) = ±{e{s)). 

We introduce the lifting r of r by 

?o = |J {(y(s), e{y(s)))} C R 2 x (R/(2ttZ)). 

s 

We define the distribution 6f (y,9) on the test function 4>(y,0) by 

<<*f >>= / ds <f){y(s),6{y{s))). 
J To 

Let us now consider a smooth evolution of the oriented closed curve T t with normal 
velocity c(y,t). For any fixed t, let us call sr t the curvilinear absissa of the curve T t and 
Sf the curvilinear abscissa of r t . Then we define the distributions 

g{y,6,t)= (jjjj+J 5 ft (y,9) and K{y,6,t) = K Vt (y) ■ g(y,6,t). 

We can now state the following result (see Theorem 2.1 and Remark 2.3 in [14] and the 
original system for (g, k) in |8J): 
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Theorem 4.1 (Transport formulation of the motion by normal velocity) 

Under the previous assumptions, the distributions g, k satisfy the compatibility equation 



r-d y g + d 9 K = (4.26) 

and the system 

g t + d y (cng) + d e ((r ■ d y c)g) + ck = _ 

K t + dyicnn) + dg((r ■ d y c)K) — (n ■ d y c)K — (r <g> r : d yy c)g = 0. \ ■ ) 

Remark that Theorem 14.11 is only known to be true for smooth evolutions of closed 
curves. When singularities appear (in general in finite time), we do not know if system 
(14271) is still true. 

4.3 Dynamics of densities of dislocation curves 

We now consider the three dimensional case with the basis (ei,e2,e 3 ). Let us now not 
only consider dislocations in the particular plane ( J/1,2/2) for y 3 = 0, but also in parallel 
planes for any y 3 = constant. To this end, let us consider distributions g(y, 9, t), n(y, 9, t) 
with the new notation y = (2/1, 2/2, 2/3) , which are again assumed to satisfy (I4.26P and 
(14.271) . To close the system we have to explain how is defined the normal velocity c(y, t) 
in the case of dislocation dynamics. To this end, we have first to compute the strain 
e(y,t) G Kgym which is solution of the system 

div a = with a = 2jie + A(trace(e))/(i 

inc e = (curl row (b x (3)) sym with 0(y,t) = J R y (27rZ) dO r(9)g{y,9,t) " ' 

where the operator inc e is obtained, taking first the curl of the column vectors of the 
matrix e, and then the curl of the row vectors of the new matrix. The curl ro ^ is the 
curl of the row vectors of the matrix, and the index ( ) sym means that we consider the 
symmetric part of the matrix. The quantity b Cg> (3 is called the Nye tensor of dislocation 
densities, where b is the Burgers vector of the dislocations under consideration. Then, 
to close the model, we define the normal velocity field by 

c{y,t) = (b®e 3 ):o-{y,t). (4.29) 

Relation (14.291) between the normal velocity and the right hand side (called the resolved 
Peach-Koehler force) is the simplest expression when the drag coefficient is isotropic and 
taken equal to 1. Again, if we neglect the short range dynamics, the complete system is 
fl4. 26ft - (j4T27l )- (14.281 )- (14. 291 ). and can be interpreted as a mean field model, describing the 
densities of dislocation curves moving in interaction. 

Finally, remark that the Groma, Balogh system <\1M is a particular case of this 
complete model, when every quantity is invariant in 2/2 and x\ — y%, X2 — 2/3, with the 
notation 9 ± (x 1 ,x 2 ,t) = g(y 1 ,0,y 3 ,±7T,t). 
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5 Concluding remarks 



The models presented in this paper deal with the dynamics of dislocation densities, and 
can be seen as a very first (and naive) attempt to derive elasto-visco-plasticity behaviour 
of crystals. All these models are essentially mean field models, and for instance are not 
able to explain relations between the stress and the plastic deformation velocity like for 
instance the typical power-law behaviour 

e p ~ ±H m 

with m large. More realistic models should also take into account the complicated short 
range dynamics, with the possible pinning of dislocations, the creation and destruction 
of dipoles of dislocations with opposite Burgers vector, and even the possibilities of junc- 
tions between dislocations, or the Frank-Read sources and the cross-slip phenomenon. 
At least, it seems that the framework of kinetic equations for the creation and annihi- 
lation of such dipoles should be a promising tool for the modelling of the short range 
dynamics of dislocations. We hope to investigate this modelling in a future work. 
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